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Data from e + e~ annihilation into hadrons collected by the JADE experiment at centre-of- 

mass energies between 14 and 44 GeV were used to study the 4-jet rate using the Durham 
\^») | algorithm as well as the first five moments of event shape observables. The data were compared 

, with NLO QCD predictions, augmented by resummed NLLA calculations for the 4-jet rate, 

f"^ ' in order to extract values of the strong coupling constant as. The preliminary results are 

I/-) \ a s {M z o) = 0.1169 ±0.0026 (4-jet rate) and as(A 7 /zo) = 0.1286 ± 0.0072 (moments) consistent 

with the world average value. For some of the higher moments systematic deficiencies of the 

QCD predictions are observed. 
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D \ 1 Introduction 

The production of hadrons in e + e annihilation allows precise tests of the gauge theory of strong 
' interactions, Quantum Chromodynamics (QCD). In this paper recent and preliminary analyses 

of JADE data using the 4-jet rate based on the Durham algorithm [1] and using the first five 
moments of event shape observables are presented [2,3]. 

The data used in our analyses were collected at centre-of-mass (cms) energies y/s = 14.0, 



22.0, 34.6, 35.0, 38.3 and 43.8 GeV between 1981 and 1986 with the JADE detector [4]. The data 
samples consist of 0(1000) events at y/s = 14.0, 22.0, 38.3 and 43.8 GeV while at y/s = 34.6 
(35.0) GeV about 14000 (21000) events are used. 

The software employed to perform the analyses includes the original JADE detector simu- 
lation, event reconstruction and event display programs, see [2, 3] for details. It is possible to 
use recent Monte Carlo event generators such as PYTHIA, HERWIG or ARIADNE to generate 
simulated events, to pass these through the JADE detector simulation and to reconstruct them 
in essentially the same way as the data. The event generators were used with parameter settings 
obtained by OPAL after adjusting to LEP 1 data. The original data are only available after a 
further step of data reduction has been performed, resulting in information about 4-vectors of 
reconstructed particles and some quantities for event selection. 

The selection of well reconstructed tracks from the tracking detectors and clusters from the 
electromagnetic calorimeter as well as the selection of hadronic events follows previous analy- 
ses. Most importantly requirements on visible energy and momentum, balance of momentum 
along the beam direction and track multiplicity suppress events from two-photon interactions, 
r production and other backgrounds [4-7]. 

The contribution of e + e~ — > bb events is subtracted from the data using simulated events. 
Then the data are corrected for the effects of detector resolution and acceptance of selection 



cuts by correction factors determined from simulated events before and after the JADE detector 
simulation. 

Experimental systematic uncertainties of the analyses include variation of the event selection 
cuts, comparison of data sets resulting from different versions of the JADE event reconstruction 
program and comparing experimental corrections derived from different event generators. 

Additional systematic uncertainties are studied for the extraction of values of the strong 
coupling constant as- Hadronisation uncertainties are evaluated by using different Monte Carlo 
generators to compute the hadronisation corrections. Theoretical systematic uncertainties are 
found by changing the renormalisation scale \x of the QCD predictions from n = yfs to fj. = y/s/2 
and fj, = 2yfs. 



2 4-Jet Rate 



Jets are reconstructed in the hadronic events using the Durham jet clustering algorithm [1]. 
Figure n (left) shows the data for the 4-jet rate corrected for experimental effects as a function 
of y cu t measured at yfs = 35 GeV. The data are compared to predictions from the event gener- 
ators PYTHIA, HERWIG and ARIADNE. We find good agreement between the data and the 
predictions within the statistical and experimental errors and conclude that the models can be 
used to correct for hadronisation effects. 




Figure 1: (left) The figure shows the corrected data for the 4-jet rate at *Js = 35 GeV compared with predictions 
from Monte Carlo models. The error bars represent statistical and experimental errors added in quadrature [2]. 
(right) The figure shows results for as from the 4-jet rate as a function of y/s. The error bars show the statistical 
and total uncertainties. The full and dash-dotted lines indicate the current world average value of as{M z o) [9]. 
The results at y/s = 34.6 and 35 GeV have been combined for clarity. Results from OPAL and ALEPH are shown 

as well [2]. 



The QCD predictions are next-to-leading-order (NLO), i.e. 0(a|) in leading-order (LO) 
with 0(ag) radiative corrections, combined with resummed next-to- leading-logarithm (NLLA) 
calculations [8]. We find good agreement between data and theory for large values of y cu t = 
0(1O~ 2 ) where mostly 3- and 4-jet events are found and the theory considering at most five 
partons is expected to be reliable. The resulting values of as are shown in figure ^ (right). 

The values for a s (M z o) from the data at -y/s = 22.0, 34.6, 35.0, 38.3 and 43.8 GeV are 
combined taking into account correlations between experimental, hadronisation and theoretical 
systematic uncertainties. The fits at y/s = 14.0 GeV have large experimental and hadronisation 



uncertainties and are therefore excluded from the average. The result is as(M z o) = 0.1169 ± 
0.0004(stat.) ± 0.0012(exp.) ± 0.0021(had.) ± 0.0007(theo.) , a s (M z o) = 0.1169 ± 0.0026 (total 
error), consistent with the current world average as(M z o) = 0.1182±0.0027 [9]. As shown in [8] 
the theoretical uncertainty rises with small y cut but is small within the fitranges. 



3 Moments of Event Shape Observables 



The first five moments of the distributions of the event shape observables 1 — T, C, B^, Byj, j/23 
and Mh are calculated according to (y n ) = J ymax y n 1/<t dcr/dy dy', where y denotes one of the 
observables, y max is the kinematically allowed upper limit of the observable and n = 1, . . . , 5. 

The calculation of perturbative QCD predictions in NLO (0(as) m LO with 0(ag) radiative 
corrections) involve a full integration over phase space. This analysis is thus complementary to 
tests of the theory using the differential distributions which are commonly only compared with 
data in restricted regions, where the theory is able to describe the data well, see e.g. [5]. 
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Figure 2: (left) The corrected data for the first five moments of the observables 1 — T, Mh and C are presented 
with error bars showing statistical and experimental error added in quadrature. The lines indicate the predictions 
of Monte Carlo models. The lower panels show the differences between data and models at y/s — 14 and 35 GeV, 
divided by the errors [3]. (right) Measurements of as(M z o) using fits to moments of six event shape observables 
are shown. The inner error bars represent statistical errors, the middle error bars include experimental errors and 
the outer error bars show the total errors. The dotted line indicates the weighted average described in the text; 
only the measurements indicated by solid symbols were used for this purpose [3]. 



Figure |2] (left) presents the data for the first five moments of 1 — T, Mh and C corrected 
for experimental effects compared with predictions by the same event generators as in sectional 
There is generally good agreement between data and model predictions; HERWIG is seen to 
describe the data somewhat less well than PYHTIA or ARIADNE. We will use the models to 
derive hadronisation corrections in order to compare the data with predictions from perturbative 
QCD. 

We fitted the QCD predictions corrected for hadronisation to the data for a given observable 
and moment n = 1, ... ,5 individually with ag(M z o) as the only free parameter. The results 
for as(M z o) are summarised in figure [21 (right). The fit to (Mh) did not converge and therefore 



no result is shown. We observe that the values of as(M z o) increase with n for the observables 
((1-T) n ), (C n ) and (B%), while for the other observables (Bfy), ((y 2 s) n ) and (Afg), n = 
2, . . . , 5, the results are fairly stable. 

We evaluated the ratio K of NLO and LO coefficients for the six observables used in our fits 
and found a clear correlation between the steeply increasing values of as(M z o) and increasing 
values of K with n for ((1 -T) n ), (C n ) and (BJ). The other observables (Bfy), ((jm)") and 
(M||), n = 2,..., 5, have fairly constant values of K and correspondingly stable results for 
as(-Mzo). We also noted that (Mh) has a large and negative value of K which is the cause that 
the fit did not converge. 

In order to find a combined value of aa,{Mip) we considered only those results for which 
the NLO term is less than half the LO term (i.e. \Kas/2ii\ < 0.5), namely (1 — T), (C), (-Bt)i 
(Byj) and ((1/23)™), n = 1,...,5 and (-Mg), n = 2,..., 5; i.e. results from 17 observables in 
total. The purpose of this requirement was to select observables with an apparently converg- 
ing perturbative prediction. Correlations between statistical, experimental, hadronisation and 
theoretical uncertainties were considered when forming the average. The result is as(M z o) = 
0.1286 ± 0.0007(stat.) ±0.0011(exp.) ± 0.0022(had.) ± 0.0068(theo.) , a s (M z o) = 0.1286 ± 0.0072 
(total error), above but still consistent with the world average value. It has been observed pre- 
viously in comparisons of distributions of event shape observables with NLO QCD predictions 
with renormalisation scale fi = \fs that fitted values of a${Mjp) tend to be large, see e.g. [10]. 

4 Summary 

We have presented preliminary results of measurements of the 4- jet rate based on the Durham 
algorithm and the first five moments of event shape observables using JADE data at y/s = 14.0 
to 43.8 GeV. The predictions of the Monte Carlo models PYTHIA, HERWIG and ARIADNE 
tuned by OPAL to LEP 1 data were found to be in reasonable agreement with the data. The 
data have also been used to extract measurements of the strong coupling constant as(M z o) with 
the results a s (M z a) = 0.1169±0.0026 (4-jet rate) and a s (M z o) = 0.1286±0.0072 (moments) in 
agreement with the world average value. The higher moments of 1 — T, C and £?t are observed 
to yield systematically larger values of as(M z o). 
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